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We investigate the communication sequences of millions of people through two different chan- 
nels and analyze the fine grained temporal structure of correlated event trains induced by single 
individuals. By focusing on correlations between the heterogeneous dynamics and the topology of 
egocentric networks we find that the bursty trains usually evolve for pairs of individuals rather than 
for the ego and his/her several neighbors thus burstiness is a property of the links rather than of 
the nodes. We compare the directional balance of calls and short messages within bursty trains to 
the average on the actual link and show that for the trains of voice calls the imbalance is signif- 
icantly enhanced, while for short messages the balance within the trains increases. These effects 
can be partly traced back to the technological constrains (for short messages) and partly to the 
human behavioral features (voice calls). We define a model that is able to reproduce the empirical 
results and may help us to understand better the mechanisms driving technology mediated human 
communication dynamics. 



INTRODUCTION 

Egocentric networks consist of a central node (the ego) 
and its immediate neighbors. They are broadly used 
in psychology and sociology as they are crucial in un- 
derstanding the interactions between an individual and 
his/her proximate social circle [TJ [2] . Dunbar discovered 
that the number of social ties of an individual is anthro- 
pologically determined to be around 150 [3H5] and that 
the ties can be classified into self-containing intimacy 
circles [6 based on the strength of one's relationships. 
These differences in the importance of the ties of an ego 
are reflected by the dynamics of social relationships. An 
individual interacts heavily only with very few people like 
close kins or closest friends and communicates less with 
most of his/her acquaintances depending on their social 
closeness and common tasks [7]. This diversity is natu- 
rally present also in the communication dynamics thus 
data collecting communication sequences of individuals 
can even be used to measure the intensity of relation- 
ships [H[T0]. Figjlja illustrates the communication pat- 
tern in an egocentric network, where the overall activity 
of the ego (green row) and activities on separated edges 
with three friends (orange rows) are presented. After 
building up the complete social network from egocentric 
subgraphs it turns out that strong links with heavy com- 
munication are usually found inside densely connected 
groups or communities as those people who are connected 
through an important relationship share many common 
neighbors. At the same time weak links more likely con- 
nect communities [8UTU]. which makes them important 
for processes evolving in social networks [TT] . 

Individual human behaviour shows heterogeneity not 
only in topology but also in dynamics. The actions of 
a person are not evenly distributed but rather clustered 



in time (e.g. see Fig{T]a). It commonly happens that 
several events are executed in bursts within a short time 
frame and such high activity periods are separated by 
long inactive ones. This kind of temporal inhomogene- 
ity can be observed in various human activity sequences 
ranging from communication to library loans or printing 
logs [T2HT5] and they are typically characterized by the 
broad distribution of inter-event times P(U e ) ~ t~ 7 . In 
the past few years different mechanisms were proposed to 
describe the origin of these heterogeneities [16 -20 , and 
the role of external influences e.g. circadian patterns have 
also been addressed 21 . 
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FIG. 1. Call activity in an egocentric network, (a) Over- 
all activity: darker the color larger the number of calls at the 
given hour. First line (green) denotes the overall in/out com- 
munication of the selected user, while the lines ei,e2 and e3 
demonstrate the communication on each active edges. (b,c,d) 
In- (red) and outgoing (blue) activity on given edges during 
selected time periods. The width of the coloured strips reflects 
the length of the corresponding calls. The selected periods are 
between 8AM and 1AM on the 5th, 16th and 29th days. 
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Even though the inter-event time distribution assigns 
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temporal inhomogeneities, it does not characterize pos- 
sible correlations in the sequences. It was observed that 
such correlations between bursty events, indeed, exist as 
indicated by trains of bursty events in Fig(l]b-d, where 
the incoming and outgoing calling activity of an individ- 
ual are shown on each of his/her links for some selected 
time periods. By appropriately defining such correlated 
trains it turned out that the distribution of the number 
E of events they contain shows a P(E) ~ E~^ scaling 
[22] . This behavior has been found in various human 
communication sequences and even in earthquake and 
neuron firing statistics, indicating same kind of univer- 
sality behind these correlated heterogeneous dynamics. 
These phenomena were interpreted as a result of memory 
effects and modeled [22] by simple multi-scale reinforce- 
ment processes similar to that describing face-to-face in- 
teractions [23] [24] . 

In spite of some related works [20] the fine grained 
structure of correlated trains in communication se- 
quences are still mainly unexplored. In this paper by 
using detailed electronic records of large number of indi- 
viduals we aim to study this problem and addressing the 
following questions: How are the event trains correlated 
with the skeleton of the backgrounding social network? 
Are they results of collective behavior or are they induced 
by individuals? Are there any correlations between the 
directions of the interactions? This paper is organized as 
follows. After introduction we define the basic quanti- 
ties and show that the correlated bursty communication 
pattern can be attributed to links rather than to nodes. 
In the subsequent section we compare the balance of in- 
and outgoing calls and SMS-s on a link to that of the 
trains on the same link. We then present a model study 
to gain deeper insight into the governing rules of the ob- 
served behavior. After drawing the conclusions we close 
the paper by the detailed description of the investigated 
datasets and used methods. 



DYNAMICS IN EGOCENTRIC NETWORKS 

In [22] it has been demonstrated that the correlations 
in a bursty time series showing a power law inter-event 
time distribution P(U e ) [T3HT5] can be detected by look- 
ing for bursty event trains and investigating the distri- 
bution P(E) of the number E of events they include. A 
train of bursty events consists of consecutive events all 
with inter-event times ti e < At inside the train, separated 
from the rest by £j e > At. For any independent P(U e ) 
(including power law ones) P(E) will decay exponentially 
while for correlated signals a power-law behavior of P(E) 



Phone calls 



P(E) ~ E~ 



(1) 



appears to indicate strong correlation of bursty events. 
In the present case for voice calls (SMS) both P(t ie ) and 
P{E) are distributed like a power-law with exponents 
7 ~ 0.7 (7 ~ 1.0) and f3 = 4.2 (f3 = 3.9) (see original solid 
lines in Figj2|. This behaviour is remarkably different 
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FIG. 2. Characteristic distributions of call and SMS 
sequences, (a) Distributions of E length of bursty peri- 
ods of outgoing events of nodes towards all the neighbors 
(solid lines) and to the direction of single neighbors (dashed 
lines). The measurement was repeated with time windows 
At = 100, 300, 600 and 3600 seconds using the original event 
sequence (original) and after the inter-event times were shuf- 
fled (random). Inset: Inter-event time distributions between 
consecutive outgoing events of the same user towards all 
neighbors (solid line) and between outgoing events on a sin- 
gle link (dashed line), (b) The same measurements repeated 
for outgoing short messages with At = 60, 120, 300 and 600 
seconds. 



from the P(E) in sequences where inter-event times of 
the whole data set are randomly shuffled. This procedure 
leaves the inter-event time distribution and the number 
of communication events between pairs of users unaltered 
but destroys all temporal correlations. Correspondingly, 
the P(E) distribution decays exponentially (random solid 
lines on Fig{2ja and b). 

If we assume burstiness to be a mere consequence of in- 
dividual human behavior, it is natural to concentrate on 
outgoing voice calls and SMS'. However, beyond study- 
ing only the timings of events, we can look at how they 
are distributed in the egocentric topology. Intuitively one 
may assume that correlated outgoing calls or SMS' of an 
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individual serve the information processing of a group, 
thus these events are directed towards several neighbours 
leading to the evolution of larger temporal motifs. The 
existence of such behaviour was demonstrated in [26] and 
also in Fig{T]c where beyond the ego two other friends 
are involved in some bursty sequences. Supposing this 
mechanism be relevant then the burstiness would be the 
property of a single node or a group of individuals. 

Surprisingly, the generic picture seems to be different. 
If we assume that the correlated events in a train are 
directed toward several neighbours, trains of events on 
single edges between two persons should have an entirely 
different, less correlated statistics of bursts. However, 
this is not the case as the P(E) distributions which were 
detected on single edges are scaling similarly and can be 
characterized by the same exponent values as in the case, 
when calls on any egocentric edges were taken into ac- 
count. This suggests that trains of events usually evolve 
on single links. Such behaviour was confirmed both for 
calls and SMS' (Figj2ja and b), where the P(E) distribu- 
tion of trains evolving between pair of individuals (origi- 
nal dashed lines) are only slightly different from the P(E) 
distributions of trains which can involve several people 
(original solid lines). This picture is also supported by 
the statistics of temporal motifs [26] , where motifs involv- 
ing two individuals are by far the most common ones. 

The same conclusion can be drawn by counting the 
number of neighbours BN whom an individual called 
in a bursty train of E events. If a user communicates 
with only one neighbour during the period then the ratio 
BN/E = 1/E or, e.g., if each calls are directed toward 
different neighbours than BN/E = 1. The distributions 
of the BN/E ratios for each E together with the aver- 
age values are presented in Figj3]a and b for calls and 
SMS, respectively. For shorter trains the distributions 
show some dispersion, however the mean values confirm 
that usually only one or two people are called in a bursty 
train. To estimate what fraction of trains of the same 
size involves one, two or more people, we calculate the 
cumulative fraction of trains involving BN number of 
neighbours. As it is depicted in Fig|3]c and d, for any 
E the vast majority of the trains of events involve only 
one neighbour (red areas) while the fractions with mul- 
tiple call receivers are considerably smaller. In addition 
Figj3]e and f indicate that even though trains consist of 
events executed with several friends, most of the calls 
are performed with only one of them. This can be con- 
cluded as the average of the maximum number of calls 
performed with the same neighbour in a train with size 
E goes with the train size as (max(CN))(E) ~ E. 
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FIG. 3. Measures of single link bursty trains, (a, b) 

The distribution and average values of the ratio BN/E (BN 
being the number of neighbours, whom an individual called 
in a bursty train) for each E train size. Trains were detected 
with At = 600 (At = 300) for calls (SMS). The pointed and 
solid lines denote the limiting case 1 and 1/E, while dashed 
line belongs to 2/E. (c, d) Cumulative fraction of the number 
of bursty neighbours (BN) for trains with size E for calls 
and SMS. (e, f) Average of the maximum number of events 
directed to the same user within a bursty period in case of 
calls and SMS. 



lated actions which are clustered into long bursty trains, 
and which trains are usually evolving on single links 
rather then connecting a larger group of people. Now 
we are at the position to ask further questions about the 
dynamics of human interactions on a dyadic level like: 
What is the direction of events in the trains? Is the com- 
munication in bursty periods balanced or dominated by 
one partner? Are there differences between voice calls 
and SMS' from this point of view? 

For the entire recording period of the dataset one 
can calculate the overall communication balance for each 
edge e as: 



max(NA, Nb) 
N a + Nb 



(2) 



MUTUAL BURSTY BEHAVIOUR 

In the previous section we pointed out similarities in 
the dynamics of communication between calls and SMS'. 
In both cases the temporal distributions of events are 
strongly heterogeneous, the sequences consist of corre- 



where Na (Nb) assigns the total number of calls from 
A to B (B to A). Hence b e can vary between 1/2 (com- 
pletely balanced) and 1 (completely imbalanced, domi- 
nated by one of the participants). Strong difference oc- 
curs between voice call and SMS dynamics if we look at 
the balance of edges and their correlation with the length 
of trains which evolve on them. To do so, we calculate 
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FIG. 4. Directed balance measures in call and SMS 
sequences, (a) Average b e edge balance values calculated 
for trains with the same size for calls (orange circles) and 
SMS (brown squares), (b) Average p(E) train balance values 
for trains of the same size in case of calls (red circles) and 
SMS (blue squares). A similar average was calculated for cor- 
responding independent event trains (yellow circles for calls 
and green squares for SMS), (c) Average p(E) and p indep (E) 
measured in call sequences where events with duration smaller 
than d> were removed, (d) The difference of the empirical 
and independent values in call sequences where events with 
duration smaller than d> were removed. 



the weighted average of b e over all trains of length E: 



bAll(E) = (b e )E trains 



(3) 



where n e (E) means the number of trains of length E 
on edge e. In Figj4ja for SMS (brown squares) b^u is 
converging to 1/2 indicating that the trains evolve on 
strongly balanced links and this effect is enhancing for 
longer trains. The balanced communication in SMS se- 
quences has been observed by Wu. et al. [20 and can 
be explained by the technologically determined way of 
communication exchange, which requires mutuality for 
SMS. However, this constrain does not apply for the mo- 
bile calls (orange points in Figj4ja) where during a call 
information can flow in both directions. Here bAii reflects 
strongly unbalanced communication and it increases to- 
wards 1 for trains with larger E. It assigns an opposite 
trend for calls compared to SMS communication as longer 
trains evolve on more unbalanced edges. 

The next question we should address is whether the 
balance or imbalance within the trains simply reflect the 
bias due to b e of the actual edge or additional effects play 
role in there. To study this we define the balance within 
a train in a similar way as for an edge: 



Pe(Em) 



max(riA, tib) 
riA+riB 



(4) 



Here p e (E m ) is the balance of the m— th train of length 
E m on edge e connecting A and B, ua (tib) denotes the 



number of events initiated by A (B) towards B (A) in 
that train; E m = tla + tib- Averaging over trains of 
the same size p(E) = (pe(E m )}E trains gives an estimate 
for the average communication balance in trains of size 
E (note that in this case different p e (E m ) values can 
evolve even for trains on the same edge e). This has to 
be compared to the case, when a train is composed from 
events selected in an independent manner from a set with 
balance b e . The latter can be calculated as 



indep 



1 E 



+ - K(i-b e y 



(5) 

where the first factor after the summation weighting the 
binomial distribution taking into account that the im- 
balance can evolve in both directions, i.e., parallel or an- 
tiparallel to the imbalance of the edge. As p L ^ ldep {E rn = 
E) depends only on b e and E, the average can be taken 
similarly to p(E) as p^P(E) = (p™** (E m )) E trains to 
get the estimate for the independent case. 

Fig. |4|b shows p(E) and p inde P(E) for both voice calls 
and SMS messages. The first apparent feature is that for 
large values of E the b A u{E) of Figjija and p inde P(E) be- 
come very similar, corresponding to the fact that for long 
trains p inde P{Ei) approaches the bias of the actual edge. 
However, the interesting effect is the difference between 
p(E) and the corresponding p znde P(E). It shows that 
trains of calls (red points) are much more unbalanced 
than one would expect from independent processes (yel- 
low circles) considering the bias of the edge. At the same 
time for SMS the contrary is true as trains (blue squares) 
are much more balanced than one would derive from ran- 
dom processes (green squares) applying the b e values of 
the corresponding edges. This demonstrates real corre- 
lation between events of the same train and suggests dif- 
ferent correlated mechanisms behind call and SMS dy- 
namics. One may suspect that the strengthening of the 
imbalance for voice calls may be a consequence of call 
trials picked up by the answering machine. To exclude 
this we repeated the calculations after removing events 
with duration shorter than <i> as depicted in Fid4lc. The 
calculated p di ^(E) = p(E) - p inde P (E) in Figjljd shows 
that the saturation value is slightly decreasing as we re- 
move longer calls, however, the differences remain signif- 
icant and saturating to 0.1 with <i> > 4 seconds. This 
duration time also serves as a timescale after which the 
answering machine calls do not play any significant role. 
For further details about the effect of short calls see Data 
and Methods. 



MODEL 

In this section we introduce a model, which is able 
to reproduce the empirical observations, namely the en- 
hanced imbalance of the bursty trains for voice calls and 
the opposite behaviour for the SMS message compared to 
independent processes. Our aim is to identify the differ- 
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ent mechanisms controlling the dynamics of communica- 
tion through these channels, integrate them into a single 
model which we can test against the empirical results. 

The emergence of correlated bursty trains in commu- 
nication sequences were interpreted earlier by memory 
processes [22] . The evolution of trains with size scaling 
as Eq{l]were explained as the result of reinforcement dy- 
namics where the probability q(n) to perform one more 
event in a train after n events have been executed, de- 
pends on n as: 

q{n) = f — -j-y J where v = /3 + 1. (6) 

Assuming this dynamics and by fitting (3 to empirical 
data, we are able to generate model trains with realistic 
size scaling. However, the question remains how we can 
introduce the mechanisms responsible for the enhance- 
ment of bursty trains. 

Let us first concentrate on voice calls. One correlation 
we observed (see Figj4ja lower panel) is that longer trains 
tend to be more unbalanced, meaning that they are more 
dominated by one of the callers. We also disclosed the 
possibility (in Figjljb) that this behaviour is an artifact 
due to answering-machine calls where the caller repeti- 
tively recalls the friend to pass an important information. 
Keeping in mind that mobile calls enable bidirectional 
information change, we assume that the observed unbal- 
anced communication in call trains reflects the difference 
in motivation between the communicating partners. If 
there is a task to solve, which is more important for one 
party, it gives motivation for him/her to repeated calls 
until the issue gets settled. 

This mechanisms can be incorporated into the rein- 
forcement process of bursty trains in the following way 
(see Fig{5]a). We simulate bursty trains, which evolve on 
a link between a pair of individuals A and B. To initiate 
a train with a probability equal to b e we randomly select 
A or B who then perform one event towards the other 
agent and at the same time we set the actual train size 
to n = 1. The decision about the next event is carried 
out in two steps. First we decide with the probability 
in Eqj6] whether to perform one more event in the train 
or initiate a new train otherwise. If the train should be 
continued a new decision is made about the direction of 
the call. The probabilities that it is initiated by A or B 
are given as: 

Ti Ti 

q*{n\(Ji) = — — or q a (n\^a 1 ) = l — (7) 

n + 1 n + 1 

where a G {A, B}. Here q a {n\<Ji) denotes the probabil- 
ity that the nth event of the actual train is performed by 
the same user who initiated the train at n = 1, while 
q a (n\^<Ji) gives the probability that the other player 
makes the call. Consequently, the longer a train evolves, 
the larger is the probability that the agent, who initiated 
the actual train will make a call to the other agent. Eqj6] 
and Eq(7| capture the coupled reasons for the evolution 
of long unbalanced trains. They reflect the reinforced 



motivation of an individual induced by the effort what 
he/she already invested in the actual series of calls to 
successfully solve a task with the other partner. 




E 



FIG. 5. Model definition and results, (a) Illustrative def- 
inition of the model. Here events are simulated between two 
people A and B. The dynamics and direction of the events 
are determined by probabilities q(n), qA(n) and gs(n) denned 
in the text, (b) P(E) distributions of empirical call trains 
(red circles) with At = 600s on edges with 0.5 < b e < 0.55 
and in corresponding model trains (black circles). The same 
functions are shown for SMS trains (blue and yellow squares 
accordingly), (c) Balance values calculated for empirical call 
(red circles) and SMS (blue squares) trains and in correspond- 
ing model processes (red and yellow symbols). Balance values 
of independent trains are also shown (green triangles) calcu- 
lated by Eq[5]with b e = 1/2. 

In case of SMS the mechanism for developing strong 
balance in bursty trains is different. There, in single 
events information can pass only one way and consec- 
utive events in a train usually have reversed direction. 
Trains are possibly conversations between partners and 
the technology constrain is responsible for the strongly 
balanced communication. To simulate this behaviour 
again we generate trains as in the previous case by using 
Eqj6j however, to select the direction of the actual event 
we define a different mechanism. Here we assume that 
the direction of an event depends only on the direction 
of the previous one as they tend to be reversed. Accord- 
ingly, conditional probabilities can be used to decide the 
direction of the actual event as: 

77/ TI 

<?a(^hcr n _i) — — , q a (n\cF n -i) = 1 — (8) 

n + 1 ri + 1 

where a G {A, B} and q a (n\^a n -i) denotes the proba- 
bility to choose the opposite direction for the nth event 
compared to the one in the n — 1th step. Accordingly 
Qa(n\(Tn-i) denotes the probability of choosing the same 
direction as for the previous event. In this way the longer 
a train evolves, the larger is the probability to revert 
the direction of consecutive events and consequently the 
more the train becomes balanced. 
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The evolution of enhanced balance/inbalance in trains 
can be checked in the best way on edges where the over- 
all communication is completely balanced and the bal- 
ance/unbalance of trains is induced only by actual be- 
havioural differences. In this way by setting b e = 1/2 
the results obtained from the model process can be com- 
pared to averages calculated for real trains that evolve on 
edges with the same b e values. We selected edges with 
overall balance values between 0.5 < b e < 0.55 and cal- 
culate the corresponding P{E) and p(E) functions for 
this limited number of 115,277,534 calls and 69,288,504 
SMS. As it is shown in Figj5]b, the size of call trains (red 
circles) detected with At = 600s and SMS trains (blue 
squares) with At = 300s are distributed broadly and 
were characterized by an exponent f3 = 4.6 and /3 = 3.5 
accordingly. Also the p(E) balance values calculated for 
the limited event sets in FigJEJc show similar behaviour 
to what was observed earlier in Figj4]b for averages cal- 
culated for all edges. This figure demonstrates that even 
if the overall communication between people is balanced 
yet the strong communication unbalance for voice calls 
and an enhanced balance for SMS trains evolves. This 
is even more conspicous when we compare these curves 
to the values obtained for independent processes (green 
triangles in Fig{5]c) using Eq{5]with b e = 1/2. 

Using the above parameters we executed model pro- 
cesses for calls and SMS with the same number of events 
and corresponding v exponents deduced from f3 accord- 
ing to Eqj6j The P(E) distributions calculated for the 
evolving model trains are fitting well to the corresponding 
empirical functions as it is depicted in FigjEJb for model 
call trains (black circles) and for model SMS sequences 
(yellow squares). It should be emphasized that there is 
no further fitting parameter in the model, nevertheless, 
the average p(E) balance values calculated for trains in 
the model processes are in very good agreement with the 
empirical observations in FigjEJc. This indicates that the 
assumed mechanisms are capturing rather accurately the 
salient features of the dynamics of directed human com- 
munication through phone calls and SMS. The only dis- 
crepancy is for the p(E) values of short SMS trains, where 
the empirical data show an even-odd effect, which is not 
present in the model, indicating that for such communi- 
cations an additional mechanism may be present enhanc- 
ing the tendency towards the balance. 



CONCLUSION 

In this paper we have studied mobile phone call and 
SMS communication sequences of millions of individuals 
and found the signature of internal correlations, which 
evolve between events and are responsible for long cor- 
related event trains. By considering egocentric networks 
we realized that these bursty trains are likely to evolve 
between pairs of individuals rather than characterizing 
the communication of an individual in a larger group. 
Moreover, after a careful analysis we have found that the 



communication in such trains is much more balanced (un- 
balanced) for SMS (voice calls) than expected from trains 
of independent events with the same balance value b e of 
the edge, where the train was formed. We have concluded 
that the backgrounding mechanisms for the evolution of 
bursty trains can be interpreted as memory processes and 
implemented by reinforcement model dynamics. We have 
shown that events in trains are determined by different 
backgrounding motivations for calls and SMS; clearly in 
the latter case they are strongly mediated by technolog- 
ical constrains. We integrated these mechanisms into a 
model to simulate call and SMS train sequences. The 
model is capable in reproducing the observed statistics 
of events in trains of bursts and create correlations to 
enhance balance (imbalance) for SMS (voice call) cas- 
cades. We believe that this work gives insight into the 
communication dynamics of individuals, dyads, and ego- 
centric networks at the level of single events. It points 
out similarities and differences of communication through 
different channels and contributes to our understanding 
of human behaviour. 



DATA AND METHODS 

In the present study we investigate the dynamics of 
human communication by analysing sequences of mobile- 
phone calls (MPC) and short messages (SMS) induced 
by a large set of individuals. The datasets were recorded 
by a single operator with 20% market share in an undis- 
closed European country. They contain 633,986,311 time 
stamped MPC (209,316,760 SMS) events recorded dur- 
ing 182 days with 1 second resolution between 6,243,322 
(4,819,993) individuals who are connected via 16,783,865 
(10,339,274) edges. In order to take into account only 
true social interactions and avoid commercial communi- 
cation, we used only actions which were executed on links 
between users, who are at least once mutually connected 
each other during the recorded period. Note that due to 
the partial coverage we cannot know the complete ego- 
centric network of each individual. This may play an 
effect on our results, however if this impact would be sig- 
nificant it would be already visible from the results of the 
limited dataset. 

To characterize temporal behaviour of an individual, 
we commonly used inter-event times defined as U e = 
ti+i — (ti +dj) where t{ and denote the starting time 
of two consecutive outgoing call events, while d{ is the 
corresponding duration. Naturally, in case of outgoing 
SMS the same definition holds but with d{ = 0. 

To filter the artificial effect of technology those con- 
secutive short messages, which were sent to the same 
neighbour with inter-event times smaller than 10 seconds 
were replaced with a single event as they were possibly 
parts of the same multipart SMS [25 j. A technology re- 
lated problem for voice calls could be that long trains 
are formed, when the caller can reach an answering ma- 
chine only instead of the callee and repeats the call sev- 
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FIG. 6. Effect of calls with short duration, (a) Cumu- 
lative distribution of duration of calls in bursty trains with 
length E. (b) P(E) distributions of periods including events 
on single edges of any direction. Distributions were calcu- 
lated for sequences where calls with duration smaller than d> 
were removed. The two straight line are denoting power-law 
functions with exponent f3 — 4.1 and 5.0. 



eral times resulting in a series of rather short calls. This 
effect is observable in Fig(6ja, where the cumulative dis- 
tribution Psr(d(E)) of call duration for calls detected in 



trains with size E is shown indicating that the length of 
calls in longer trains are shorter. In order to check the 
effect of such calls on the statistics the distributions were 
calculated after calls with duration d > d> were removed 
(Figj6|. Altough the exponent changed from /3 = 4.1 to 
5.0 the power-law scaling of P(E) remained indicating 
strong residual temporal correlations. 
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